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Guru P Guruswamy 
Informatics General Corporation 
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SUMMARY 

The development and application of transonic smal 1 -di sturbance codes at NASA Ames Research Center for 
coiUputing two-dimensional flows, using the code ATRAN2, and for computing three-dimensional flows, using 
the code ATRAN3S, are described Calculated and experimental results are compared for unsteady flows about 
airfoils and vnnqs, including several of the cases from the AGARD Standard Aeroelastic Configurations In 
two dimensions, the results include AGARD priority cases for the NACA 64A006, NACA 64A010, NACA 0012, and 
MBB-A3 airfoils In three dimensions, the results include flows about the F-5 wing, a typical transport 
wing, and the AGARD rectangular wing Viscous corrections are included in some calculations, including 
those for the AGARD rectangular wing For several cases, the aerodynamic and aeroelastic calculations are 
compared with experimental results 


SYMBOLS 

b semichord of wing 

c ful 1 chord of wing 

Cj^ sectional lift coefficient due to plunging mode 

Cj sectional lift coefficient due to pitching mode 

sectional moment coefficient due to plunging mode 

sectional moment coefficient due to pitching mode 

k reduced frequency based on chord 

? semispan of the wing 

M,M^ free-stream Mach number 
U flutter speed 

X sectional distance, measured in semichords, from the elastic axis to the mass center 

a 

u wing-section-to-air-mass density ratio 

(j) disturbance velocity potential 

bending natural frequency 
reference frequency 
cj pitching natural frequency 


1 INTRODUCTION 

In the last decade, there have been extensive developments in computational, unsteady transonic aero- 
dynamics (Refs 1 and 2) This growth in computational methods was in response to a need by engineers for 
computer codes with which to study fundamental aerodynamic and aeroelastic problems in the critical tran- 
sonic regime For example, because the transonic flight regime provides the most efficient aircraft cruise 
performance (Refs 3 and 4), transonic flow fields constitute one of the most intensely studied problems in 
fluid dynamics Most large commercial aircraft cruise in the transonic regime However, the current total 
cost (Ref 5) of developing a modern transport aircraft is so large that it puts most of the resources, as 
well as the future reputation, of a company at risk Computational fluid dynamics (CFD) provides a new 
tool that in combination with the use of test facilities such as wind tunnels, provides aerodynamic design 
information with which to reduce this risk at the earliest possible phase of aircraft development 

One of the major tasks in developing a new transport that involves unsteady transonic flow is the 
flutter analysis of the supercritical wings Experiments have shown (Ref 6) that dips in the flutter 
boundaries for wings occur at transonic Mach numbers and that such dips are especially severe for super- 
critical wings This occurrence of instabilities at lower dynamic pressures in the transonic regime is 
attributable to the motion of the shock waves that are present on the wings (Ref 7) Proper modeling of 
the physics of such moving shock waves requires that the CFD methods solve nonlinear partial differential 
equations for regions of mixed subsonic and supersonic flow At the present time, the most advanced codes 
use the smal 1 -di sturbance , transonic potential equation, these codes are being used for generic research 
in aeroelasticity Currently, codes are being developed that use the more exact full-potential equation 
In comparison, for steady flows, practical applications use full-potential codes, and there is extensive 
development of Euler codes 
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The successful development and application of the two-dimensional code LTRAN2 (Ref 8), which solves 
the unsteady, transonic, smal 1 -di sturbance potential equation using a time-integration method, and the 
availability of faster computers with larger memories made possible the development and use of the three- 
dimensional code XTRAN3S (Ref 9) In this code, the alternating direction implicit (ADI), finite- 
difference scheme (Ref 8) in LTRAN2 was extended to three-dimensional flow over wings In this conserva- 
tive method, shock waves move in a time-accurate manner The XTRAN3S code can also perform static and 
dynamic aeroelastic computations by simultaneously integrating the aerodynamic and structural equations 
of motion This code is being used for both aerodynamic and aeroelastic applications (Refs 10-13) 

Several efforts are in progress to develop two-dimensional and three-dimensional codes based on a 
full-potential theory (Refs 14-20) The reason for going to the more exact ful 1 -potent i al formulation 
from the approximating, smal 1 -di sturbance one is to remove the accuracy limitations in the latter The 
removal of these limitations may be necessary to adequately model the flow over supercritical wings, and 
also to handle large perturbations in the flow However, by using the ful 1 -potenti al formulation, the 
major advantage of simplicity (imposing boundary conditions inherent in the smal I -disturbance theory) is 
lost Modeling a complex geometry (such as a wing, body, and nacelle) is fairly easy in smal 1-di sturbance 
theory, for the fu 1 1 -potential theory, however, the resultant computational meshes are quite complex 
This feature is one reason why full-potential codes are not as well developed as the smal 1-di sturbance 
potential codes 

This report describes the development and application of transonic smal 1 -di sturbance codes at NASA 
Ames Research Center for computing two-dimensional flows, using the code ATRAN2, and three-dimensional 
flows, using the code ATRAN3S Comparisons of calculated and experimental results are presented tor 
unsteady flows about airfoils and wings, including several of the cases from the AGARD Standard Aeroelastic 
Configurations In two dimensions, the results include AGARD priority cases for the NACA 64A006, 

NACA 64A010, NACA 0012, and MBB-A3 airfoils In three dimensions, the results include flows about the 
F-5 wing, a typical transport wing, and the AGARD rectangular wing Viscous corrections will be included 
in some calculations, including those for the AGARD rectangular wing For several cases, the aerodynamic 
and aeroelastic calculations are compared with experimental results Brief descriptions are also given 
of recent improvements of the algorithms used in these codes, included are the addition of high-frequency 
terms and more accurate differencing in ATRAN2 , and grid improvements in ATRAN3S for f i ghter-ai rcraft-type 
wi ngs 


2 ALGORITHM DEVELOPMENTS AND FLOW-FIELD CALCULATIONS TWO DIMENSIONS 

In this section we shall describe algorithm and code improvements to the NASA Ames code LTRAN2 and 
then present comparisons of computed results and experimental data for the AGARD airfoils (Ref 21) The 
LTRAN2 code was introduced 7 years ago (Ref 8) Since then, its use has become routine (Refs 1 and 2) 
in generating unsteady, transonic air loads for oscillating airfoils, particularly for use in applications 
to aeroelasticity (Kefs 22 and 23j During those 7 years, there have been many improvements in the accu- 
racy, capability, stability, and efficiency of LTRAN2 Improvements in accuracy include the addition of 
high-frequency terms in the boundary conditions (Ref 24) and the governing equation (Ref 25), the inclu- 
sion of viscous effects (Refs 26-28), and the use of second-order-accurate spatial differencing (Ref 29) 
Wind-tunnel wall modeling (Ref 30) and supersonic free-stream (Ref 31) capabilities were added The 
stability of the code was improved by using monotone algorithms (Ref 32), and, finally, the efficiency of 
the code was improved by using nonreflecting far-field boundary conditions (Ref 33) 

The improvement in stability with the monotone algorithm is shown in Fig 1 (Ref 32) The older 
method can become unstable at sonic expansion points and hence require a much smaller time-step for stable 
calculations (12 times smaller in this case) 

The improvement in accuracy with second-order spatial differencing is shown in Figs 2 and 3 (Ref 29) 
for flow with type-B shock-wave motion Because of the large excursion in the shock wave, which the 
second-order method describes more accurately (Fig 2), the unsteady pressures (Fig 3) are more accurately 
described over a large portion of the airfoil by the second-order method 

Figure 4 (Ref 25) shows the improvement in accuracy with the inclusion of high-frequency terms 
This modification did not require the addition of a third level of computer storage, as did previous 
methods, which can be an advantage especially when the method is extended to three-dimensional algorithms 
With the inclusion of these high-frequency terms in LTRAN2, the code is now valid at all frequencies 
This version of the NASA Ames code is designated ATRAN2 code 

To account for viscous effects, two procedures are available in ATRAN2, a viscous-ramp method and a 
lag-entrainment method (Refs 27 and 28) For inviscid calculations, ATRAN2 solves the unsteady, transonic, 
smal 1 -di sturbance potential equation (Refs 25 and 32) 

+ 2kM^+ t = [(1 - M^) - (y + +4. (1) 

The viscous corrections are implemented by modifying the inviscid airfoil tangency conditions, as well as 
the downstream wake conditions in the lag-entrainment method The viscous-ramp procedure is a phenomeno- 
logical method in which a prion determined shape-changes simulating the viscous displacement effects are 
incorporated into the inviscid procedure Here the interaction between the shock and boundary layer is 
modeled by placing a wedge-nosed ramp at the base of the shock to obtain the reduced shock pressure rise 
The lag-entrainment equations are based on the boundary-layer assumption that the normal extent of the 
VISCOUS region is small when compared with airfoil or wake thickness By integrating the governing partial 
differential equations in the normal direction and suitably modeling the requisite relationships, a set of 
three first-order ordinary differential equations, lag-entrainment equations," is obtained References 27, 
28, 34, and 35 provide details on these methods and their use in ATRAN2 
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2 1 NACA 64A010 Airfoil 

See Ref 34 for details on the procedures for generating the computed results Unsteady inviscid, 
wedge, and lag-entrainment solutions were obtained by forcing the airfoil to undergo sinusoidal motion and 
then integrating the flow equations in time, starting from the corresponding steady-state solutions A 
Fourier analysis of all three solutions was performed for pressure distributions for the third cycle of 
oscillation, and the real and imaginary parts of the first Fourier component were extracted These results 
were plotted after defining the pressure as 

Cp = j , (Re SI n t - Im cos .t) 

where is the amplitude of oscillation Unsteady aerodynamic coefficients were computed, based on the 

third cycle, for selected values of reduced frequencies Using these coefficients, flutter speeds were 
computed by the U-g method, following the procedure described in Ref 34 Two degrees of freedom, 
pitching and plunging, were assumed for the aeroelastic model 

In Refs 36 and 37, steady and unsteady aerodynamic results for the NACA 64A010 airfoil, obtained in 
the Ames 11- by 11-Foot Transonic Wind Tunnel, are given for various flow conditions In this study, 
results are computed and compared with wind-tunnel results at a Reynolds number of 12 6 x lo® for Mach 0 796 
and a mean angle of attack of -0 21 The airfoil configuration is taken from the AGARD report (Ref 21) 

Figure 6 shows the lower-surface, steady-pressure distributions obtained by computations and the 
experiment The viscous-ramp and the lag-entrainment methods compare better in predicting the shock loca- 
tion than does the inviscid method Unsteady aerodynamic results were ODtained for the pitching motion of 
the airfoil at the five reduced frequencies - 0 05, 0 102, 0 202, 0 302, and 0 404 - that were considered in 
the experiment The inviscid and viscous-ramp calculations used 720 time-steps per cycle, the lag- 
entrainment method required 8000 In Fig 6, the plots of the magnitude and the corresponding phase angle 

of the lift coefficient versus the reduced frequency are shown, and Fig 7 shows the magnitude and phase 
angle of the moment coefficient about the leading edge as a function of the reduced frequency 

For k = 0 101 and 0 404, Fourier analyses of all three solutions were performed for lower-surface 
pressure distributions for the third cycle of oscillation, and the real and imaginary parts of the first 
Fourier component were extracted A comparison of these results with experiment is shown in Figs 8 and 9 
In general, the viscous solutions, particularly the lag-entrai nment solution, are closer to the experimental 
results than is the inviscid solution 

Using the U-g method, flutter speeds were computed based on the unsteady aerodynamic coefficients 
obtained by the computations and the experiment They were computed by assuming the elastic axis at the 
quarter chord, a = -0 5 , the mass center at the midchord, x^ = 0 5, and the ratio of pjunge-to-pitch 
natural frequencies " 0 ^ the nondimensionalized flutter speed 0, and the corre- 

sponding reduced frequency k, are shown versus the airfoil-to-air-mass density ratio, y Flutter speeds 
obtained by the lag-entrai nment solution are closer to those of the experiment than are the other two solu- 
tions In general, all three theoretical methods compare fairly well with the experiment in the reduced 
frequency, it is also observed that viscous effects have a tendency to reduce the flutter speed 

Figure 11 is a time-history of the computed lower-surface pressures for the inviscid flow during the 
third cycle of motion Notice the motion of the shock wave and the change in shock strength To simulate 
the shock motion correctly, the governing equation was solved in conservative form Also notice the oscil- 
lations in pressures at the leading edge, which are due to the use of airfoil coordinates from the experi- 
mental model (Ref 36) 

2 2 The MBB-A3 Airfoil 

For this airfoil (Ref 34) results were computed and compared with those of experiment (Ref 38) at a 
Reynolds number of 6 x 10® for a Mach number of 0 765 and a mean angle of attack of 1 5° (the design condi- 
tions for the airfoil) As was done in Ref 38, the Mach number and the angle of attack considered for 
computations were 0 7557 and 1 30°, respectively, thus matching the flow conditions in the wind tunnel 
The airfoil configuration is taken from the AGARD report (Ref 21) 

Figure 12 shows the comparison of the steady-state pressure distributions for the upper and the lower 
surfaces The lag-entrai nment solution compares well with the experiment The wedge solution is closer to 
the experiment than the inviscid solution Steady lifts obtained by inviscid, viscous ramp, lag-entrai nment, 
and the experiment are 0 6667, 0 6433, 0 5377, and 0 5190, respectively 

Unsteady aerodynamic results (Ref 34) were obtained by pitching the airfoil about the midchord with 
an amplitude of 1 0° for four reduced frequencies, 0 05, 0 10, 0 15, and 0 20 The unsteady inviscid, 
viscous-ramp, and lag-entrainment solutions were obtained starting from the corresponding steady-state 
solutions obtained earlier For all the cases considered, three cycles were required during which the 
transients disappeared and a periodic solution was obtained The inviscid and viscous-ramp calculations 
used 720 time-steps per cycle, and the lag-entrainment used 8000 

For k = 0 1, a Fourier analysis of all three solutions was performed for the upper-surface pressure 
distributions for the third cycle of oscillation, and the real and imaginary parts of the first Fourier 
component were extracted A comparison of these results is shown in Fig 13 Using the U-g method, 
flutter speeds (Ref 34) were computed based on the unsteady aerodynamic coefficients obtained by the three 
methods The viscous methods increased the flutter speeds, but the reduced frequencies were essentially 
the same as those for the inviscid method 

2 3 NACA 006 Airfoil 

Three AGARD cases (cases 6, 8, and 10 in Ref 21) were computed for the NACA 006 airfoil and compared 
with experimental data (Ref 37) The computations were made using the experimental parameter given in 
Ref 37 
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For case 6, the \/alues used were = 0 853, the dynamic flap amplitude -Sq ~ 1, k = 0 12, and the 

mean flap angle ^ 0 16 number of time-steps per cycle was 720 for all the remaining two- 

dimensional results Only the inviscid and viscous-ramp methods v/ere used for the remaining results 
because of the extremely small time-step size required for the 1 ag-entrai nment method 

Figure 14 shows the comparison of pressures for the lower surface for steady flow, and Fig 15 shows 

the unsteady comparisons for the lower surface For all the remaining cases, the unsteady computational 

results were taken from the third cycle of motion Figure 16 shows a time-history of the lower-surface 
pressures for the inviscid results for the third cycle of flap motion Notice the motion of the shock 
wave and the change in its strength It essentially disappears during a portion of the cycle, that is, 
type B shock-wave motion This case demonstrates the need for a time-accurate formulation of the computed 
results Finally, notice the weak ridge and valley caused at the three-quarter chord location owing to 
the flap motion 

For case 8 (Ref 21), = 0 877, k = 0 118, 6 q ^ 1 13, and = 0 15 Figure 17 shows the compari- 

sons of steady pressures for the lower surface, and Fig 18 shows the unsteady flow comparisons The 
steady calculations predict too strong a shock wave, and this overprediction is carried over to the unsteady 

results, which show the effects of the shock motion As noted in Ref 37, these differences can be par- 

tially ascribed to wi nd-tunnel -wal 1 interference In previous calculations (Ref 8) on this airfoil with 
an oscillating flap, similar differences were found between calculations and experiment, for the cases of 
type A, B, and C shock-wave motion For those cases, calculations that included wi nd-tunnel -wal 1 simula- 
tions (Ref 30} also indicated that the differences might be due to wind-tunnel-wall interference 

Figure 19 shows a time-history of the inviscid results Notice the motion of the shock wave and the 
changes at the flap hinge location Figure 20 shows the variation in the lift coefficient. Cl, during the 
third cycle of motion In comparison with harmonic approximations, zero through the first harmonic and 
zero through the fourth harmonic. Fig 20 shows that the variation is described by the first harmonic 
Figures 21 and 22 show the variations in lower shock position and strength during the third cycle, where 
the shock strength is taken to be the difference in pressure coefficients across the shock wave The 
shock-position variation is essentially simply harmonic, whereas the shock-strength variation contains 
higher harmonics This variation in shock position requires a method that allows the shock to move over 
the airfoil 

For case 10 (Ref 21), = 0 879, k = 0 468, 5 q = 1 08, and 6^ = 0 01 Figure 23 shows the com- 
parisons of steady pressures and Fig 24 shows the unsteady flow comparisons for the lower surface As in 

case 8, the calculations overpredict the shock strength in the steady flow, and the effects are carried 
over into the unsteady flow comparisons Notice the effects of increasing the frequency in comparing 
Figs 18 and 24 At the higher frequency, the shock motion is less and the real part of the pressure 
coefficients has a dip upstream of the shock wave Figure 25 shows a time-history of the inviscid results 

For this relatively high frequency, there is little shock-wave motion, as can be seen by comparing 

Figs 19 and 25 

2 4 NACA 0012 Airfoil 

For the NACA 0012 airfoil, most of the cases (Ref 21) required large mean angles of attack or large 

changes in the angle of attack that are typical for helicopter applications The only case that appears 

to be within the capability of the small -disturbance formulation is case 5 As we did for the previous 
airfoil, we have chosen to use the experimental parameters (Ref 37) for computing this case to compare 
with experimental data For case 5, = 0 755, k = 0 1628, the mean angle of attack - 0 016°, and 

the dynamic pitch angle riQ = 2 51 ° The airfoil is pitching about the quarter-chord and that is where 

the moment coefficients are computed 

Figure 26 shows the comparison of steady pressures for the upper surface Figure 27 shows comparisons 
at the eight times, during the cycle of motion, for which the experimental data (Ref 37) were given 
Figure 28 shows a time-history of the inviscid results For this high dynamic pitch angle, Fig 28 shows 
the creation, large motion, and disappearance of a strong shock wave Figures 29-33 show the variations, 
during the third cycle, of the lift coefficient Cl, of the quarter-chord moment coefficient Cm, of the 
drag coefficient Cd, of the upper shock position, and of the upper shock strength Figure 29 shows that 
Cl is described adequately by the first harmonic Figure 30 shows that Cm requires higher harmonics 

for an adequate description Figure 31 shows that C^ is adequately described by the harmonics through 

the second, recall that for the drag, its fundamental harmonic is twice the first harmonic of the motion 
Figure 32 displays the shock-position variation and the portion of the cycle where the shock disappears 
Finally, Fig 33 shows the strong variation in shock strength, over the portion of the cycle where the 
shock disappears, there are no marks for actual shock strength 

No results are shown for the NLR 7301 airfoil Calculations showed poor agreement with experiment 
Because of the thickness (16 5^) and leading-edge bluntness of this airfoil, it is probably outside the 
capability of the smal 1 -di sturbance formulation 


3 ALGORITHM DEVELOPMENTS AND FLOW-FIELD CALCULATIONS THREE DIMENSIONS 

In this section, we shall describe the progress to date that has been made in a program at Ames 
Research Center for developing and applying an efficient computer code for calculating unsteady transonic 
flows in three dimensions We have chosen to develop further and to apply the XTRAN3S code, which was 
developed by Borland et al (Refs 39-41) under contract to the United States Air Force 

The choice of XTRAN3S was based on several considerations First, XTRAN3S was a direct extension of 
LTRAN2 to three dimensions Hence, the successful results produced by LTRAN2 led us to expect that 
XTRAN3S would be an efficient and accurate tool for studying three-dimensional flows Like LTRAN2, the 
governing equation that is solved by XTRAN3S is the unsteady, transonic, small -disturbance equation in con- 
servation form The use of the small-disturbance equation has the advantage of allowing the use of simple 
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computational grids for the unsteady flow about complex geometrical configurations, such as a wing, body, 
and canard, with, for example, moving control surfaces Since the algorithm relies on shock-capturing to 
account for the location and movement of shock waves, the equation is solved in conservative form, other- 
wise, erroneous shock locations and speeds will occur 

The development of the Ames version of XTRAN3S, which is called ATRAN3S, proceeded as described below, 
where each of four developments is described in a separate subsection that also includes computational 
results The first effort was an aerodynamic and aeroelastic study (Refs 10 and 11) to compare XTRAN3S 
calculations with experimental data The purpose of that effort was to validate results obtainable from 
the code, as well as to determine the proper choice of numerical parameters, such as time-step size, in 
order to obtain accurate results In that first study, the simple case of a rectangular wing was examined 

The next step (Ref 12) was to apply the code to a transport-type wing, characteri zed by a high-aspect 
ratio, large-taper ratio, and small sweep angle A flutter-analysis comparison was made with experimental 
data More recently (Refs 42 and 43), a study was made of the severe case of a fighter-type (F-5) wing, 
characterized by a low-aspect ratio, small-taper ratio, and large sweep angle Finally, a study (Ref 44) 
was made of the effects of viscosity and modes on the transonic aerodynamic and aeroelastic characteristics 
of wings The AGARD rectangular wing, as well as a typical transport wing, were examined The latter two 
studies required that algorithms be developed For the fighter wing, a new computational grid was devel- 
oped in order to stabilize the calculations In the study of viscous effects, the manner of interacting 
the viscous and inviscid calculations was modified to eliminate numerical fluctuations in the pressure 
profiles at the shock locations 

3 1 Aerodynamic and Aeroelastic Study Rectangular Wing 

In this first study, XTRAN3S calculations for a rectangular wing were compared with experimental data 
The purpose was to validate XTRAN3S 

Aerodynamic equations of motion - Many forms of the smal 1 -di sturbance equations have been developed 
for computing the transonic flow field about wings (Refs 26 and 45) In this analysis, the modified 
unsteady, three-dimensional, transonic smal 1 -disturbance equation is used 

where is the disturbance velocity potential, A = M^, B = 2M^, E = (1 - M^) , F = -(l/2)(y + 1)M^, 

G = (1/2)(y - 3)M^, and H = -(y - 1 )M^ 

For the results shown in this section, the low-frequency form of this equation was solved by setting 

A to zero and using corresponding boundary conditions. The corresponding computer code, LTRAN3, was the 

earlier version of XTRAN3S and was the only version available at the time of this study (Ref 11) The 
code IS based on a time-marching, finite-difference scheme following the first-order-accurate alternating 
direction implicit (ADI) algorithm A detailed description of the procedure is given in Ref 39 For 
details on the computational grid and the convergence criteria for the steady calculations, see Refs 10 
and 11 

Unsteady aerodynamic pressures were computed by forcing the wing to undergo a sinusoidal modal motion 
and integrating the aerodynamic equation of motion in time The modal motion assumed was the same as that 
simulated in the experiments For all the cases studied here, it was found that about three cycles of 
motion with 360 time-steps per cycle were sufficient to obtain a periodic aerodynamic response Periodicity 
was tested by comparing the responses of the second and third cycles The magnitudes and phase angles of 
the unsteady pressure jumps and corresponding force coefficients were computed using the third cycle 

Aeroelastic equations of motion - Generalized coordinates (Refs 10 and 11) were used in deriving the 
aeroelastic equations of motion In this analysis, two generalized coordinates h(t,y) and a(t,y), which 
correspond to bending displacement and torsional rotation of the elastic axis of the wing, respectively, 
were chosen as representative of the fluttering wing The generalized coordinates h(t,y) and a(t,y) can 
be expressed as 

h(y,t) = h(t)f(y) , a(y,t) = 5(t)e(y) 

where h(t) and a(t) are unknown functions of time, and f(y) and e(y) are assumed semirigid modes For 
further details on the aeroelastic equations of motion, see Refs 10 and 11 

The aeroelastic parameters and sign conventions for a typical section of the wing are shown in Fig 34 
It is assumed that the wing is rigid in the chordwise direction, and that the amplitudes of oscillation are 
small It is also assumed that the principle of superposition of air loads is valid, even in the presence 
of shocks The validity of this assumption has been shown for two-dimensional cases both by experiment 
(Ref 36) and theory (Ref 23), provided the shock wave does not introduce separation and that the shock 
motion is small 

Flutter solution procedure - A procedure based on the U-g method (Ref 46) was used to determine the 
transonic flutter boundaries Unsteady aerodynamic coefficients required in this work were the generalized 
lift and moment coefficients owing to modal motions corresponding to a pure bending mode f(y) and a pure 
torsional mode e(y) For the U-g method, unsteady aerodynamic coefficients, Cj,a» ®nd Cfna are 
required as a function of reduced frequency for each mode In this analysis, the coefficients were com- 
puted at three reduced frex^uencies , and they were interpolated by a Lagrange interpolation scheme Since 
a low-frequency assumption was used in LTRAN3, the reduced frequencies considered were less than 0 4 (based 
on full chord) 

Results Comparison of aerodynamic pressure coefficients - In Ref 47, experimental investigations 
were conducted on an unswept rectangular wing The wing had an aspect ratio of 3 with a 5%-thick biconvex 
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airfoil section Both steady and unsteady pressures were measured Unsteady pressures were measured while 
the wing was oscillating in its first bending mode 

In Ref 11, steady and unsteady pressures from LTRAN3 were compared with the experimental data at 
Mach numbers of 0 7 and 0 9 for four span stations located at 0 %, 50'“ , 70 , and 90"^ semi span For unsteady 
computations, the same bending mode that was measured in the experiment was also simulated in the code 
These results were also compared with correspond! ng data obtained in Ref 47 in which linear aerodynamic 
theory, based on the kernel function method, was used In the following results, the magnitude of the 
unsteady pressure jump is scaled by the induced angle of attack corresponding to the amplitude of the tip 
displacement, and the phase angle is defined as positive if the pressure leads the bending displacement 
The magnitude and phase angle from LTRAN3 correspond to the first fundamental harmonic in a Fourier series 
decomposition of the pressure time-history 

In Fig 35, steady-pressure curves are compared between experiment and LTRAN3 at M = 0 9 at four 
spanwise stations The two sets of curves compare fairly well Except at the root section, pressure 
coefficients and shock locations obtained by LTRAN3 are in close agreement with those obtained from the 
expenmeni; The discrepancies at the root section can be attributed mainly to the boundary layer on the 
wall, which was not considered in LTRAN3 

In Fig 36, magnitudes |ACp| and phase angles of the unsteady pressure jump obtained by LTRAN3, 
experiment, and kernel-function method are plotted at the 70‘a semi span station for M = 0 9 and k^ = 0 26 
In general, the two sets of curves obtained by LTRAN3 and experiment compare fairly well Peaks in pres- 

sure jumps occur at almost the same locations for both LTRAN3 and experiment As expected, linear theory 
results do not compare favorably, either with LTRAN3 or experiment, in the vicinity of shock 

Results Comparison of flutter boundaries - Reference 48 reports an experimental investigation of the 
transonic flutter characteristics of unswept rectangular wings with aspect ratios of 5, and with circular- 
arc sections at various thickness ratios and Mach numbers See Refs 11 and 48 for details of the experi- 
mental models 

In Refs 10 and 11 results obtained from LTRAN3 and from experiments were compared for four cases 
(1) 6 % thick at M = 0 715, (2) 6% thick at M = 0 851, (3) 6% thick at M = 0 913, and (4) 4% thick at 
M = 0 904 These cases were selected so that Mach numbers ranging from a no-shock case to a strong-shock 
case were included Based on studies using LTRAN3, the Mach numbers considered in the experiment did not 
include a moderate-shock case for the 6%-thick model Thus, a case from the 4%-thick model was selected 

Based on the unsteady aerodynamic coefficients obtained from LTRAN3, flutter boundaries were computed 
by the U-g method In Fig 37, results from LTRAN3 are plotted as a curve of flutter speed and corre- 
sponding reduced frequency versus wing-air-mass-density ratio for a 6%-thick model at M = 0 715 The 
corresponding curve obtained by the linear code NASTRAN is given in the same figure The experimental 

results available for a wing-air-mass-density ratio of 36 72 and a reduced frequency of 0 232 are also 
shown 


Similar results were obtained for the other three cases Results for the four cases are given in 
Table 1 From this table it can be seen that the flutter speeds obtained by LTRAN3 are greater than those 
obtained by experiment On the other hand, the reduced frequencies obtained by LTRAN3 are lower than 
those obtained by the experiment, except for case 1 With the increase in Mach number, both LTRAN3 and 
the experiment show an increase in flutter speed, a decrease in reduced frequency, and an increase in 
density ratio Comparisons are better at lower Mach numbers Differences are quite significant at 

M = 0 904 and M = 0 913 These differences can be mainly attributed to the discrepancies between the 
aerodynamics of LTRAN3 and the experiment rather than to the flutter modeling 

3 2 Flutter Analysis of a Transport Wing 

The next step in evaluating the code XTRAN3S was to apply it to a flutter analysis of a transport-type 
wing (Ref 12) A search of the literature (Ref 12) available at that time revealed an appropriate set of 
experimental data with which to make comparisons in a Japanese report (Ref 49) The model configuration 
is shown in Fig 38 It has an aspect ratio of 8, a taper ratio of 0 4, a quarter-chord sweep of 20°, and 
a NACA 65A012 airfoil section See Ref 12 for a description of the model's structural and mass data 

Steady and unsteady flows were calculated in the region of the transonic dip in the flutter boundary 
in the experimental data at Mach numbers of 0 75, 0 8, 0 825, and 0 85 Figure 39 shows the steady-flow 
pressure coefficients at M = 0 85 The flow is supersonic over a large portion of the wing and terminates 
in a shock at approximately the 65% chord line There were no aerodynamic data in Ref 49 with which to 
compare, although in Ref 12 comparisons were made with a standard code for steady flows and showed good 
agreement 

Unsteady aerodynamic coefficients were calculated for the first three uncoupled vibration modes They 
were the first bending, second bending, and first torsion modes for a swept-wing beam model of the struc- 
ture For each mode and Mach number, the wing was oscillated sinusoidally for three cycles at reduced fre- 
quencies that bracketed the expected flutter-reduced frequency Approximately 1200 time-steps per cycle 
were used for these calculations Figure 40 shows a comparison with linear theory for a reduced frequency 
k =0 121 (near flutter) and a Mach number of 0 85 for the first bending mode There are substantial dif- 
ferences because of the transonic effects The larger leading-edge moment coefficients from XTRAN3S indi- 
cate a rearward shift of the aerodynamic center caused by the lift contribution of the shock-wave motion 

Using these unsteady aerodynamic data, a flutter analysis was performed, using the U-g method 
(Ref 46) Figure 41 shows comparisons of flutter boundaries obtained from XTRAN3S, a doublet lattice 
analysis, and the experimental data The XTRAN3S results agree better with experiment than do the linear 
theory results Results for XTRAN3S are shown at Mach numbers of 0 75, 0 80, and 0 825 No flutter was 
predicted at Mach 0 85 
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3 3 F-5 Wing Analysis 

Analysis - For nonrectangular wings, such as the F-5 wing, XTRAN3S used a sheared coordinate system, 
in which the shearing occurs in the plane of the wing (Fig 42), where 


£(x,y) 


X - x^gCy) 
XiE(y) - ’ 


n(y) = y 


After applying the transformation, £q (2) is transformed to 


c(z) = 2 
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The finite-difference computations are done in the f,n,t coordinate system, as indicated in Fig 42 

This conventional shearing transformation, as given by Eq (3), is simple and adequate for wings with 
small sweep and large taper ratio (Refs 9-13) However, for wings with low aspect ratio, high sweep, 
and small taper ratio, such as fighter wings, use of this shearing transformation produces unstable calcu- 
lations Since the transformation given by Eq (3) is only a function of the local chord, computational 
flow regions in the physical domain depend on the wing planform For delta-type wings, Eq (3) yields 
highly skewed flow regions and thus large discontinuous values for the metrics, Cy, near the upstream- 
and downstream-flow boundaries The metric Cy appears as a coefficient of the cross derivatives in the 
governing equation, as given by Eq (4) Also, the far-field boundaries in the physical domain are not 
aligned with the flow directions These combined factors can make the flow computations unstable, as was 
the case for the F-5, AV-8B, and DAST ARW-2 wings (Ref 13) 

The flow domain generated by Eq (3) will be illustrated for the F-5 wing Figure 43 shows the plan- 
form of the F-5 wing which has an aspect ratio of 2 98, a taper ratio of 0 31, and a leading-edge sweep 
angle of 31 92° The computational flow region that was obtained by using the shearing transformation, 
given by Eq (3), is shown in Fig 44 From this figure, the large skewness of the grid lines, which 
causes large gradients for the metrics near the flow boundaries, can be observed In fact, the scale in 
Fig 44 is stretched 10 to 1 in the y (or vertical) direction So the actual skewness is 10 times worse 
than shown Also the kinks (i e , discontinuities) in the upstream and downstream boundary grid lines are 
10 times worse than shown It is at the downstream boundary at the location of the kink in the coordinate 
line (i e , at the wing tip span station) that numerical instabilities originate in the calculations for 
the F-5 wing 

A modified shearing transformation was developed (Refs 42 and 43) that eliminated these numerical 
instabilities Like Eq (3), the modified transformation maps the physical wing into a rectangle, as indi- 
cated in Fig 42 Away from the wing, the modified transformation has the following characteristics 
(1) far-field boundaries are independent of the wing planform and aligned with respect to the free-stream 
direction, (2) smooth first and second derivatives occur for values of the metric quantities, particularly 
near boundaries, and (3) grid lines are clustered near the leading and trailing edges Figure 45 shows 
the computational flow region for the F-5 wing that was obtained by the modified shearing transformation 
Detailed comparisons of the metric quantity ^y for the two methods are given in Ref 43 The new trans- 
formation, along with other improvements in XTRAN3S, has been incorporated into the Ames version of this 
code, which is called ATRAN3S 

Resul ts - In Ref 43, using the new transformation technique, steady and unsteady aerodynamic results 
were computed for the F-5 wing at M = 0 80, 0 90, and 0 95 and were compared with experimental measure- 
ments (Ref 50) 

Figure 46 shows plots of steady pressure results obtained at M = 0 9 for the mid-semispan station 
by the two transformation methods and experiment In spite of using a time-step size of 0 001, results 
from the conventional shearing transformation eventually diverged Also, those results were highly inaccu- 
rate after 4000 time-steps, as illustrated in Fig 46 The computations diverged sometime after 4000 time- 
steps and before 6000 time-steps With the modified shearing transformation, a converged solution was 
obtained by using 2000 time-steps of size 0 01 In Fig 46, it can be seen that the latter method compares 
well with the experiment Plots of steady pressure distributions for four span stations obtained by the 
modified transformation and the experiment are given in Fig 47 Comparisons are generally good at all 
span stations 

Figure 48 shows the modal motion used in the NLR experiment The wing is pitching about an axis 
located at the 50% root chord, and the pitching axis is normal to the wing root Figure 49 shows plots of 
the real and imaginary values of the upper-surface pressures at four span stations obtained by the modified 
shearing transformation and the NLR experiments at M = 0 9 These results were obtained for the wing 
oscillating at a frequency of 40 Hz 

The same modal motion used in the NLR experiment was simulated in the code Results from the code 
were obtained by forcing the wing to undergo a sinusoidal modal motion for three cycles with 1200 time-steps 
per cycle, during which time the transients disappeared and a periodic response was obtained Because of 
the instabilities encountered during the steady calculations, no attempt was made to use the conventional 
shearing transformation 
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3 4 Effects of Viscosity and Modes 

Anal ysis- In the final study (Ref 44), ATRAN3S was used to examine the effects of viscous corrections 
on the aerodynamics about wings and the effects of including higher structural modes in addition to the 
fundamental bending and torsion modes in transonic aeroelastic analyses Two wings were studied, a rec- 
tangular wing, with a NACA 64A010 airfoil section, and a swept wing, with an MBB-A3 supercri tical airfoil 
section Viscous effects on both wings were analyzed by employing the viscous-wedge and lag-entrainment 

methods Aeroelastic analyses were performed and the effects of including the first two bending and tor- 

sion modes in the analysis were determined Results from this work show that the inclusion of viscous 
effects increases the flutter speed for the two wings studied For the rectangular wing, the fundamental 
modes were sufficient to determine the flutter speed, but the second torsion mode was required for an 
accurate aeroelastic analysis of the swept wing 

Viscosity can play an important role in both the aerodynamic and aeroelastic characteri sties of wings, 
as indicated in Ref 34, where the two-dimensional code, LTRAN2 (viscous) was used to study airfoils In 
transonic flow, viscous effects can alter the shock location and strength These changes in the pressure 
distribution along the chord will also influence the aeroelastic characteristics of the wing The ATRAN3S 
code was modified to include the effects of flow viscosity (Ref 41) Two models are used to account for 
the viscous effects in ATRAN3S The viscous wedge is an empirical method and the lag-entrai nment method 
IS based on a set of integral boundary-layer equations These methods are the same as those first incor- 
porated into LTRAN2 (viscous) (Refs 28 and 41) The same viscous equations used in two dimensions are 
applied stripwise along the span to yield a three-dimensional correction 

This viscous-wedge method is two dimensional and is applied to the wing at each spanwise coordinate 
Details of the computations can be found in Ref 44 A two-dimensional shock profile was used to deter- 
mine the placement of the wedge in this study, as opposed to the three-dimensional shock profile discussed 
in Ref 41 This change was made to help reduce oscillations in the pressures obtained when using the 
viscous methods 

Solutions obtained using the original viscous corrections contained high-frequency oscillations in 
the unsteady air loads An investigation into this problem determined that the shock sweep angles computed 
for the viscous-wedge corrections did not change smoothly across the span, but instead jittered in some 
span regions, as well as with time-step, at some span locations Also, it was determined that the use of 
these shock sweep angles to determine the shock location was not consistent with the method of locating 
the shock wave for use in the type-dependent differencing that is used in the inviscid calculations To 
correct this inconsistency, the shock location at each span station was determined by only using flow-field 
values along that span station This new method of determining the shock location was then used to place 
the viscous wedge on that span section As shown in Ref 44, large fluctuations in the pressure profile 
at the shock location were eliminated when the new method was used 

For some calculations, the wedge model alone is sufficient for accurate determination of the unsteady 
viscous pressure distributions For other calculations, however, a more exact analysis is required to 
model the viscous effects To do this, the wedge correction is used along with a set of integral boundary- 
layer equations, the lag-entrainment equations, to perform viscous calculations downstream of the shock 
wave, including the downstream wake In this method it is assumed that viscous regions are small relative 
to the wing or wake thickness Details of the two methods are given in Refs 28, 34, and 41 

In this study, the aeroelastic equations of motion for a cantilever wing were formulated using an 
assumed mode method (Refs 46 and 51) The superposition of the bending and torsion natural modes was 
used to describe the deflection shape of the wing at flutter It was assumed that the amplitude of oscil- 
lation was small, this was done so that the transonic aerodynamics from each mode of motion could be super- 
posed (Ref 23) For details on the aeroelastic equations of motion, see Ref 44 

A three-dimensional U-g type flutter analysis was used to obtain the flutter boundaries In this 
study, the first two bending and torsion modes were chosen to describe the motion of the low-aspect-ratio, 
rectangular, centi lever wing Previous studies have shown that the fundamental modes are sufficient to 
determine the flutter speed of a rectangular wing in subsonic flow For comparison, the first two bending 
and torsion modes were also used for the swept wing 

The generation of the aerodynamic data was performed separately from the aeroelastic analysis First, 
ATRAN3S was used to generate the steady and unsteady aerodynamic data by prescribing the wing motion 
These data were then input into the aeroelastic equations of motion to solve for the flutter speed The 
unsteady solutions were obtained by forcing the wing to oscillate sinusoidally as the aerodynamic equation 
of motion was integrated in time To obtain a steady solution, the amplitude of motion was set to zero 
The unsteady viscous calculations were restarted from corresponding steady-state solutions For unsteady 
calculations, three cycles of motion were sufficient for the transients to die out and for a periodic solu- 
tion to be obtained For most inviscid calculations, 360 time-steps per cycle were used Viscous calcula- 
tions typically required 720 time-steps per cycle for the wedge calculations and 2160 time-steps per cycle 
for the lag-entrainment calculations Unsteady aerodynamic coefficients, based on the third cycle, were 
used to compute flutter boundaries for inviscid and viscous calculations For the rectangular wing, 
unsteady calculations were made for four modes at four different Mach numbers and from three to six reduced 
frequencies for each Mach number Also the inviscid and the two viscous methods were considered Only 
limited results were computed for the swept wing 

Results Rectangular wing with NACA 64AQ1Q airfoil - One of the two wings considered in this study 
was the AGARD rectangular wing with an aspect ratio of 4 and a NACA 64A010 airfoil cross section Aero- 
dynamic and flutter results were computed for the rectangular wing at M = 0 8 and at zero angle of 
attack, which is one of the AGARD suggested flow conditions (Ref 52) In addition, calculations for the 

same wing were made at M = 0 70, M = 0 85, and M = 0 90 to determine the effects of Mach number on 
three-dimensional aeroelastic calculations First, steady-state pressure distributions were obtained for 
the inviscid, wedge, and lag-entrainment methods The steady wedge calculations were restarted from the 
steady inviscid solution, and the steady lag-entrainment calculations were restarted from the steady wedge 
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solution Figure 50 shows the steady pressure profiles for the inviscid flow at M = 0 85 Figure 51 
shows the steady pressure distribution for the three methods at the 60S semi span station and at three 
Mach numbers The code diverged when the lag-entrainment method was used at M = 0 90, so no results were 
obtained for those conditions The forward displacement of the shock and reduction in strength are char- 
acteristic of VISCOUS flow These trends are consistent with those shown in two dimensions for the same 
airfoil It can be seen that the lag-entrainment method shows no deviation from the wedge for this par- 
ticular wing 

The first and second bending and first and second torsion natural modes were selected to accurately 
model the motion of the wing to obtain unsteady aerodynamic coefficients Calculations were made at 
reduced frequencies, based on full chord, from 01 to 0 6 The maximum deflections were 3 44“ for the 

torsion modes and 10, of the chord for the bending modes The wing was constrained to pitch about the 

midchord Figure 52 shows the magnitude of the unsteady local lift coefficient and the corresponding 
phase angle versus span for the rectangular wing at three Mach numbers Results are shown for inviscid 
and for both the wedge and lag-entrainment viscous correction methods for the first torsion mode at a 
reduced frequency of 0 4 As expected from two-dimensional studies (Ref 34), the lift coefficients 
obtained by the viscous methods are lower than those calculated using the inviscid method, however, the 
magnitude of the lift using the wedge correction shows little change over the inviscid method The phase 
angles show a significant change for both viscous methods 

As mentioned earlier, the flutter boundaries were determined using a three-dimensional U-g method 
Of primary interest in this study were the aeroelastic characteristics of the rectangular wing related to 
the effects of the structural modes and flow viscosity Four combinations of the first two pure bending 
and torsion modes were considered First, only the first bending and first torsion modes were included 
Then the second modes of each type were added individually Finally, all four modes, namely, first bend- 
ing, second bending, first torsion, and second torsion, were all included in the solution For the modal 
analysis, the wing was configured so that the elastic axis and center of mass were located at 45% and 60% 
of the chord, respectively, aft of the leading edge All other aeroelastic parameters were set to repre- 
sent a typical wing For the rectangular wing, the first bending and torsion modes were sufficient to 
determine the flutter speed (Ref 44) The flutter frequency in all cases is very close to the first 

bending frequency, which explains the dominance of the first bending mode Figure 53 shows the deflection 

of the wing at flutter It closely resembles the first bending mode Many of the aeroelastic parameters 
were varied to see if any effect of the higher modes could be obtained, little to no effect was seen 
Only when very unrealistic values were used was any significant effect of modes noticed At all Mach 
numbers, the first bending and torsion modes are all that are required for an aeroelastic analysis of a 
smal 1 -aspect-ratio, unswept, rectangular wing in subsonic or transonic flow 

The effect of flow viscosity is an important factor to consider when computing the flutter boundaries 

of a wing in transonic flow Viscosity can have a noticeable effect on the aerodynamic forces acting on 

the wing and, hence, on the flutter boundaries Viscous calculations were made at three Mach numbers, 

M = 0 80, M = 0 85, and M = 0 90 For each Mach number, aerodynamic coefficients were calculated over a 
range of reduced frequencies from 01 to 0 6 and for the same four structural modes discussed earlier 
The wedge and lag-entrainment viscous methods were used At M = 0 85, the flow was transonic and the 
shock was well defined The lift coefficients calculated with viscous corrections included were slightly 
below the inviscid calculations This is the trend that was anticipated because it corresponds to trends 

shown for two-dimensional airfoils This small reduction in force, coupled with significant changes in the 

phase angles, contributes to the higher flutter speeds shown for the viscous methods in Fig 54 (Ref 44) 

A well-defined transonic dip can be seen in Fig 54, with the lowest flutter speed occurring between 
M = 0 85 and M = 0 90 

Inclusion of the viscous corrections resulted in about a 10% increase in flutter speed The increase 
in flutter speed is a direct result of the viscous corrections causing a reduction in the magnitude of the 
lift coefficient A higher flutter speed does not always follow from a reduction in magnitude of the lift 
because the change in phase angle must also be considered when one is trying to anticipate flutter trends 
The wedge correction required little more computational effort than the inviscid method, but yielded an 

increase in flutter speed, as is depicted in Fig 54 The lag-entrainment method on the other hand, 

required 6 to 8 times the computing time to obtain a converged, stabilized solution The results obtained 
using the lag-entrainment method were about the same as the wedge correction provided It is noted that 
the wedge calculations were performed at 720 time-steps per cycle, and the lag-entrainment calculations 
were performed at 2160 time-steps per cycle Both of these are substantially less than were used in the 
earlier two-dimensional calculations reported in Ref 34 

Results Swept wing with MBB-A3 airfoil - This swept wing with an MBB-A3 airfoil section was chosen 
to represent a typical transport-type wing It has a full-span aspect ratio of 8, a taper ratio of 0 4, a 
leading-edge sweep angle of 25°, and an MBB-A3 supercritical airfoil section Because of the complex geom- 
etry, the computing time required per case was twice that required for similar calculations on the rectangu- 
lar wing Aerodynamic and flutter calculations were made inviscidly, as well as by using the wedge and 
lag-entrainment viscous corrections at a Mach number of M = 0 85 As was done for the rectangular wing, 
steady-state data were computed, and then the unsteady runs were restarted from the corresponding steady- 
state data files 

The first bending, first torsion, second bending, and second torsion modes were incorporated into the 
flutter analysis of the swept wing (Fig 55) Again, the deflections of the wing were computed for each 
mode at 507 wing locations to describe the deflection shapes in the transonic code It can be seen that 
the modes, despite the sweep and taper, are very similar to the corresponding modes for a cantilever beam 
To obtain a converged invisci ) solution, 720 time-steps per cycle were required, 1440 and 2880 time-steps 
per cycle were required for the viscous wedge and lag-entrainment methods, respectively Calculations were 
made at k=04, 05, 06, and 0 7 for all three methods 

Figure 56 shows the inviscid steady-pressure distribution across the swept wing for M = 0 85 The 
steady pressure across the chord, at the 60% semi span location, is shown in Fig 57 for all three methods 
and for both the upper and lower surfaces The viscous wedge shows a shock which has moved forward of the 
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inviscid shock and is much weaker The lag-entrainment correction causes the shock to move still farther 
forward For this wing and airfoil section, the lag-entrainment method does show a significant change 
from the viscous wedge method 

Flutter calculations were performed on this wing at M = 0 85, the results are shown in Table 2 
Here again, the primary modes dominate the flutter solution For this wing, however, the second torsion 
mode does cause a noticeable change in the flutter speed When the second bending mode was added to the 
primary modes, only an insignificant change in the flutter speed was noted, but when the second torsion 
mode was added, the flutter speed decreased by about 5% It was expected that the torsion modes would 
play an important role in this analysis because of the sweep of the wing The first two torsion modes can 
be seen along with the first bending mode in Fig 58 


4 CONCLUDING REMARKS 

Two computer codes, ATRAN2 and ATRAN3S, were used to compute unsteady transonic flows about airfoils 
and wings, and the results were compared with experimental data Both codes solve the unsteady, transonic 
small -di sturbance equation in conservative form, and viscous corrections are made through the use of 
either the viscous wedge or 1 ag-entrai nment methods The ATRAN2 code computes two-dimensional flows about 
airfoils, and ATRAN3 computes three-dimensional flows about wings 

The comparisons with experimental data showed good agreement, in general, and included several of the 
test cases from the AGARD Standard Aeroelastic Configurations For some cases, aeroelastic calculations 
were made, and the transonic dip in the flutter boundary was demonstrated For the transport-type wing, 
experimental flutter data were available, and the calculated results were in good agreement with them 
For some airfoil cases, including the NACA 0012 airfoil, during the airfoil motion there were large varia- 
tions in the location and strength of the shock waves These cases demonstrated the need for codes to 
solve the governing equations in conservation form in order to accurately simulate the unsteady flow 

Several improvements in the algorithms were described, including improvements for high-frequency 
accuracy, numerical stability, and second-order spatial differencing accuracy for ATRAN2, the grids for 
fighter-type wings and vi scous-i nvi scid algorithm interactions for stability for ATRAN3S were also improved 

These calculations demonstrate that ATRAN2 and ATRAN3S are ready for use by aerodynamicists and aero- 
elasticians in practical applications involving unsteady transonic flows for which the limitations of small 
disturbance theory are valid 
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Table 1 Comparisons of flutter speed and corresponding reduced frequency between 

LTRAN3 and experiment 


CASE 

MC*CH NO 

THICKNESS 

RATIO 

% 

DENSITY 

RATIO 

m/7r/>b^ 

REDUCED FREQUENCY 

tJc/U 

LTRAN3 EXPT 

FLUTTER SPEED 

U/ba;^ 

LTRAN3 EXPT 

1 

0 715 

6 

36 72 

0 250 

0 232 

4 30 

3 83 

2 

0 851 

6 

58 72 

0 120 

0 162 

5 60 

4 55 

3 

0913 

6 

74 65 

0 045 

0 122 

8 80 

4 94 

4 

0 904 

4 

75 17 

0 085 

0 138 

6 60 

3 70 


Table 2 Flutter speeds for various modal 
combinations swept wing, M = 0 85 


MODES 

INVISCID 

WEDGE 

LAG 

ENTRAINMENT 

FIRST BENDING 
FIRST TORSION 

2 1237 

2 8507 

2 2066 

FIRST AND SECOND 
BENDING 
FIRST TORSION 

2 1240 

2 8504 

2 2057 

FIRST BENDING 
FIRST AND SECOND 

2 0666 

2 9511 

2 1856 

TORSION 




FIRST AND SECOND 
BENDING 

FIRST AND SECOND 

2 0670 

2 9502 

2 1850 

TORSION 






Fig 1 Comparison of algorithms using the mono- 
tone Godunov switch LTRAN2-MG and nonmonotone 
switch LTRAN2, plots of upper-surface pressure 
coefficients of a NACA 64A010 airfoil (experi- 
mental model ordinates) in pitching motion 
= 0 80 


Fig 2 Comparison of first- and second-order methods 
for unsteady flow about a NACA 64A006 airfoil with 
oscillating trailmg-edge flap, type-B shock-wave 
motion, instantaneous pressure-coefficient plots for 
upper surface kt = 213® 
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Cp = Cp srn(kt - o) 
= 1 k^t = 3" 

= 0 854 k = 0 358 


2ND ORDER 1ST ORDER 


ClI = 3 30 
Ol = 50° 

C|V| = 1 05 
0M = -I7r 


iCli =3 31 

ol ^ 

C(y] = 1 08 
0|V1 = -180° 




• PRESENT RESULTS (REF 25) 
□ BASIC LTRAN2 (REF 8) 
A LTRAN2HI (REF 24) 

LINEAR THEORY 

PITCHING MOMENT 




Fig 3 Unsteady upper-surface pressure-coefficient Fig 4 Real and imaginary components of pitching- 
plots (amplitude and phase) of the first harmonic moment coefficients for an airfoil oscillating in 
component for the case in Fig 2 pitch = 0 7, a = sin u.t 


NACA 64A010 AIRFOIL 
= 0 796 Qq =-0 21 R = 12 6 X 10® 

INVISCID 

VISCOUS RAMP 

LAG ENTRAINMENT 

. EXPERIMENT 



Fig 5 Comparison of lower-surface steady- 
pressure distributions, theories and experiment 
NACA 64A010 airfoil , M = 0 796 


NACA 64A010 AIRFOIL 
M^ = 0 796 Qq =-0 21 R = 12 6 X 10® 



40 r 


■D 30 


I 20h 


< 10 





I 1 l_ I I 

0 12 3 4 

REDUCED FREQUENCY 

Fig 6 Comparison of plots of lift coefficient 
versus reduced frequency, theories and experiment 
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NACA 64A010 AIRFOIL 
= 0 796 Oq =-0 21 R = 12 6 X 10® 



REDUCED FREQUENCY 

Fig 7 Comparison of plots of moment coefficient 
about leading edge versus reduced frequency, 
theories and experiment 


NACA 64A010 AIRFOIL 

= 0 796 Qq = -0 21 k = 0101 R = 12 6 / 10® 

INVISCID 

VISCOUS RAMP 

LAG ENTRAINMENT 


• EXPERIMENT 



012345678910 

x/c 

Fig 8 Real and imaginary parts of first Fourier 
component of unsteady lower-surface pressure coeffi- 
cients k = 0 101 
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NACA 64A010 AIRFOIL 

= 0.796, a-Q = -0.21”, k = 0.404, R = 12.6 X 10® 


INVISCID 

VISCOUS RAMP 

LAG-ENTRAINMENT 




Fig. 9. Real and Imaginary parts of first Fourier 
component of unsteady lower-surface pressure 
coefficients: k = 0.404. 


NACA 64A010 AIRFOIL 
M^ = 0.796, Q'o = -0.21”, R - 12.6 X 10® 




Fig. 10. Comparison of flutter speeds, theories and 
experiment. 


MBB-A3 AIRFOIL 


LOWER SURFACE PRESSURES CYCLE 3 
NACA 64A010 AIRFOIL 

M^ = 0.796, = -0.21°, k = 0.404, R = 12.6 X 10® 



Moo (THEORY) = 0.7557 
Moo (EXPERIMENT) = 0.765 
ofo (THEORY) = 1.3° 

(EXPERIMENT) - 1.5° 


INVISCID 

VISCOUS RAMP 

LAG-ENTRAINMENT 

® EXPERIMENT 



Fig, 11. Time-history of lower- surface pressure 
coefficients for third cycle of motion: 

NACA 64A010 airfoil, k 0.404. 


Fig. 12. Comparison of upper and lower surface 
steady-pressure distributions, theories and 
experiment. 
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MBB-A3 AIRFOIL 

- 0.7557, - 1.3°, k = 0.1, R = 6 X 10® 

INVISCID 

VISCOUS RAMP 

LAG-ENTRAINMENT 



Fig. 13. Real and imaginary parts of first 
Fourier component of unsteady upper-surface 
pressure coefficients: k = 0.1. 


Fig. 14. Comparison of lower-surface steady-pressure 
distributions, theories and experiment: NACA 64A006, 
M = 0.853. 


NACA 64A006, LOWER SURFACE 
Moo = 0.853, k - 0.12 


INVISCID 

VISCOUS (WEDGE) 




x/c 


LOWER SURFACE PRESSURES CYCLE 3 

NACA 64A006 AIRFOIL 
Moo = 0.853, k = 0.12 



Fig. 15. Real and imaginary parts of first 
Fourier component of unsteady lower-surface 
py'essure coefficients: k = 0.12. 


Fig. 16. Time-history of lower-surface pressure 
coefficients for third cycle of motion: NACA 64A006 

airfoil , k = 0.12. 
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NACA 64A006, LOWER SURFACE 
Moo = 0-877, k = 0.118 


INVISCID 

VISCOUS (WEDGE) 

o EXPERIMENT 



NACA 64A006, LOWER SURFACE 
Moo = 0.877, k = 0.118 

— INVISCID 

VISCOUS (WEDGE) 


o EXPERIMENT 



Fig. 17. Comparison of lower-surface steady 
pressure distributions, theories and experiment: 
NACA 64A006, M = 0.877. 


-5.0 



x/c 

Fig. 18. Real and imaginary parts of first Fourier 
component of unsteady lower-surface pressure coeffi- 
cients: NACA 64A006, k = 0.118. 


LOWER SURFACE PRESSURES CYCLE 3 
NACA 64A006 AIRFOIL 
Moo = 0.877, k = 0.118 



NACA 64A006 AIRFOIL 
Moo = 0.877, k = 0.118 
CYCLE 3 

HARMONIC 

0 THROUGH 4th 

— - 0 PLUS 1st 



TIME, deg 


Fig. 19. Time-history of lower-surface pressure Fig. 20. Lift-coefficient variation during third 
coefficients for the third cycle of motion: cycle: NACA 64A006 airfoil, k = 0.118, = 0.877. 

NACA 64A006 airfoil, k = 0.118, M«, = 0.877. 
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LOWER SHOCK POSITION vs TIME 
NACA 64A006 AIRFOIL 
= 0 877, k = 0118 
CYCLE 3 


HARMONIC 

0 THROUGH 4th 

0 PLUS 1st 



) 1 ^ ^ ^ i ^ ^ • ' 

720 760 800 840 880 920 960 1000 1040 1080 

TIME, deg 


LOWER SHOCK STRENGTH vs TIME 
CYCLE 3 


HARMONIC 

0 THROUGH 4th 

--- 0 PLUS 1st 
^ ACTUAL 


ujl 
o *“ 
<o 

LL Z 

QC UJ 

3 cc 

</5 h- 
to 


M8 



TIME, deg 


Fig 21 Lower-shock-position variation during Fig 22 Lower-shock-strength variation during third 

third cycle NACA 64A006 airfoil, k = 0 118, cycle NACA 64A006 airfoil, k = 0 118, = 0 877 

= 0 877 


NACA 64A006, LOWER SURFACE 
M^ = 0 879 k = 0 468 


NACA 64A006 LOWER SURFACE 
Moo = 0879, k = 0468 


INVISCID 

VISCOUS (WEDGE) 



INVISCID 

VISCOUS (WEDGE) 



x/c 


Fig 23 Comparison of lower-surface steady pres- 
sure distributions, theories and experiment 
NACA 64A006 airfoil, M = 0 879 



x/c 

Fig 24 Real and imaginary parts of the first 
Fourier component of unsteady lower-surface pressure 
coefficients NACA 64A006, k = 0,468 
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NACA 0012, UPPER SURFACE 
= 0.755, = 0.016, k - 0,1628, = 2.51 



Fig. 25. Time-history of lower- surface pressure 
coefficients for the third cycle of motion: 

NACA 64A006 airfoil, k = 0.468, = 0.879, 


Fig. 26. Upper-surface steady pressure distributions: 
NACA 0012, M = 0.755. 
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UPPER SURFACE PRESSURES CYCLE 3 
NACA 0012 AIRFOIL 

Moo " 0.755, k = 0.1628, (Xq = 2.51, = 0.016 




Fig. 28. Time-history of upper-surface pressure 
coefficients for the third cycle of motion: 

NACA 0012, - 0.755, k = 0.1628. 


IMACA 0012 AIRFOIL 

Moo = 0.755, = 2.51, k = 0.1628, = 0.016 

CYCLE 3 

HARMONIC 



TIME, deg 


Fig. 30. Quarter-chord moment-coefficient varia- 
tion during the third cycle: NACA 0012, 

= 0.755, k = 0.1628. 


NACA 0012 AIRFOIL 

Moo = 0.755, = 2.51, k = 0.1628, = 0.016 

CYCLE 3 



TIME, deg 

Fig. 29. Lift-coefficient variation during third 
cycle: NACA 0012, = 0.755, k = 0.1628. 


NACA 0012 AIRFOIL 

Moo = 0-755, = 2.51, k « 0.1628, <x^ = 0.016 

CYCLE 3 

HARMONIC 

0 THROUGH 4th 

---0 PLUS 2nd 



TIME, deg 

Fig. 31, Drag coefficient during third cycle: 
NACA 0012, = 0.755, k = 0.1628. 
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UPPER SHOCK POSITION vs TIME 
NACA 0012 AIRFOIL 

M^ = 0 755, tto = 2 51, k = 0 1628, ° 0^® 

CYCLE 3 


z 

o 

p 

< 

o 

o 


u 

o 


< 


cc 

D 

c/5 


U1 

o. 

o. 

3 


HARMONIC 



TIME, deg 


UPPER SHOCK STRENGTH vs TIME 
NACA 0012 AIRFOIL 

Moo = 0755, Oq = 251, k = 01628, = 0016 

CYCLE 3 

HARMONIC 



TIME, deg 


Fig 32 Upper-surface shock-position variation 
NACA 0012, M„ = 0 755, k = 0 1628 


Fig 33 Upper-surface shock-strength variation 
NACA 0012, = 0 755, k - 0 1628 


5% THICK CIRCULAR ARC WING 
ASPECT RATIO = 30 
M = 09 


FLOW 




Fig 34 Definition of aeroelastic parameters for 
a cantilever wing 


LTRAN3 



x/c 

Fig 35 Comparison of steady pressure coefficients 
LTRAN3 and experiment 




OVI 'dlAII 






FIRST BENDING MODE 

DOUBLET LATTICE 

XTRAN3S 

LIFT COEFFICIENTS LEADING EDGE MOMENT COEFFICIENTS 



Fig 40 Modal aerodynamics for the Japanese transport wing M„ = 0 85, k = 0 121 
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Fig 41 XTRAN3S flutter boundary for the Japanese Fig. 42. Sheared coordinate system used for 

transport wing swept and tapered wings 
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F5 WING PLANFORiVl 
AR = 2.98, TR = 0.31 



FAIRING (0.0250 m) 

F1g. 43. Dimensions of the F-5 wing. 


CONVENTIONAL SHEARING TRANSFORMATION 


y 



X 


Fig. 44. Physical grid (64 x 20) in x-y plane from 
the conventional shearing transformation. 


MODIFIED SHEARING TRANSFORMATION 




Fig. 45. Physical grid (64 x 20) in x-y plane 
by the modified shearing transformation. 


Fig. 46. Effect of transformation on steady pressures 
at M = 0.90. 




"ri3: 
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F5 WING 
AR = 2 98 
TR =0 31 
= 0 90 
a =0° 
t =40 Hz 


UNSTEADY Cp 

REAL 

IMAGINARY 

□ REAL 
A IMAGINARY 


ATRAN 3S 
NLR EXPERIMENT 


REAL 


IMAGINARY 



M = 085 

O!o = 0 



_i » > 

14 16 18 


Fig 49 Comparison of unsteady pressures 
theory and experiment M = 0 90 


Fig 50 Rectangular wing steady pressure distribu- 
tions M = 0 85 
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CHORD LENGTH 


INVISCID 

• WEDGE 

LAG ENTRAINMENT 

FIRST TORSION MODE 
60% SEMISPAN STATION 



Fig 51 Rectangular wing steady pressure shown Fig 52 Magnitude of unsteady lift and corresponding 

for three Mach numbers and the inviscid and vis- phase angle versus span for the inviscid and viscous 

cous methods at 60% semi span station and a = 0 methods 


M = 0 85 



AO 


Fig 53 The shape of the rectangular wing at 
flutter IS composed of the first bending and 
torsion modes 


INVISCID 

WEDGE 



MACH NUMBER 

Fig 54 Effect of Mach number and viscosity on 
flutter speed of the rectangular wing 
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Fig 58 Shape of the swept wing at flutter is composed of the first bending and the first and second 

torsion modes 
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